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Introduction 


Abstract 


The dwarf sheoak Allocasuarina humilis (Otto & A.Dietr.) 
LASJohnson (Casuarinaceae) is a dioecious wind-pollinated 
woody shrub endemic to the biodiverse Southwest Australian 
Floristic Region (SWAFR) of Western Australia. It occurs 
throughout most of the SWAFR, where it is found in vegetation 
communities ranging from heath to forest (Western Australian 
Herbarium 1998-). Plants grow to 2 m high, resprout from 
a lignotuber following fire and are thought to be very long- 
lived. Plants are diploid (2 n = 20; Barlow 1959), flower from 
May to November (Western Australian Herbarium 1998-) and 
infructescences contain multiple glabrous samaras 5-6 mm long 
with a short truncate wing (Wilson and Johnson 1989) that may 
facilitate some wind dispersal. Seed dispersal via emus has also 
been noted (Calvino-Cancela etal. 2006). 


Eight polymorphic microsatellite markers 
were developed for Allocasuarina humilis (Otto 
& A.Dietr.) LASJohnson, a common woody 
shrub endemic to south-west Australia that 
is a common component of the understorey 
across a range of vegetation types. The markers 
were isolated using 454 shotgun sequencing 
and tested on 48 individuals from two natural 
populations. We found six to 19 alleles per 
locus (mean 12.1 j and observed heterozygosity 
ranged from 0.208 to 0.956 (mean 0.710). The 
markers will be used to assess population 
structure and diversity in A. humilis, with the 
aim of developing appropriate seed sourcing 
strategies for ecological restoration. 
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Allocasuarina humilis is a common understorey 
component of ecological restoration programs for 
areas disturbed by logging, mining and roadworks and 
is also used in community-based restoration projects. 
Seed sourcing has been mostly conducted in areas 
local to restoration sites, following the principles of local 
provenance and the guidelines of the previous Western 
Australian Forest Management Plan (Conservation 
Commission of Western Australia 2004). However, with 
a fully outcrossing breeding system and the potentially 
widespread dispersal of pollen via wind, local sourcing 
of seed may not be the most appropriate strategy for A 
humilis. 

Molecular analysis can reveal patterns of gene flow, 
genetic diversity and population connectivity to better 
inform the seed sourcing decision-making process. 
Restoration outcomes are greatly improved when 
information on a species' genetic characteristics is 
incorporated into the planning phase of restorations 
(Godefroid et al. 2011). Here, we developed eight high 
resolution microsatellite markers for A humilis with 
which to evaluate its genetic diversity and population 
structure, to assist in developing appropriate seed 
sourcing strategies for ecological restoration. 

Materials and methods 

Genomic DNA was extracted from freeze-dried 
leaf material from one A humilis plant using the 
Doyle and Doyle (1987) CTAB method with 1% w/v 
polyvinylpyrrolidone (MW 40000) added to the 
extraction buffer. Shot-gun sequencing was conducted 
at the Australian Genomic Research Facility (Brisbane, 
Australia) using 4 pg DNA on a GS-FLX Titanium 
machine (Roche, 454 Life Sciences) following Gardner 
et al. (2011). We used QDD 1.3 (Meglecz et al. 2010) 
with default parameters to screen raw sequences for £ 
eight di- to hexa-base repeats. Primers were designed 
using Primer3 within QDD. Of the 539 unique loci, 10 
were chosen for polymerase chain reaction (PCR) trials 
following recommendations in Megleczetal. (2010) and 
Gardner et al. (2011). Loci were preferred that had pure 
microsatellite motifs £8 repeats in length that were >20 
bp from primers. 

Universal Ml3 tails (5- CACGACGTTGTAAAACGAC) 
were appended to each forward primer at the 5' end, 
and the complement sequence was modified with a 


5' fluorescent label (FAM, VIC, NED or PET). Singleplex 
polymerase chain reaction (PCR) amplifications were 
performed in an 8 pL reaction mix containing 20 ng 
template DNA (extracted from freeze-dried leaf material 
using the method detailed above), 50 mM KCI, 20 mM 
Tris-HCI (pH 8.4), 0.2 mM of each dNTP, 2.0-2.75 mM 
MgCL, 0.15 pM fluorescent labelled Ml3 primer, 0.03 pM 
forward primer, 0.15 pM reverse primer, 0.05 pL Taq DNA 
polymerase. Most reactions benefited from addition 
of 0.5 M betaine to the PCR mix (Table 1). Loci were 
amplified using one of two programs in an Eppendorf 
gradient thermocycler: (1) 95°C for 4 min; 15 cycles of 30 
s at 94°C, 30 s at 65 to 50°C with a stepdown of 1 °C per 
cycle, and 72°C for 45 s; 30 cycles of 15 s at 94°C, 15 s at 
50°C, and 45 s at 72°C; 8 min at 72°C; (2) as for Program 
1, but the initial amplification step involved 25 cycles 
with a stepdown of 0.6°C per cycle and an extension 
phase of 80 s at 72°C. Products were visualised on 8% 
polyacrylamide gels stained with ethidium bromide, 
and PCR reactions were further optimised for MgCI 2 
and/or betaine (as per Table 1). PCR products were 
visualised on an Applied Biosystems 3730 capillary 
sequencer and scored manually using GeneMapper 4.0 
(Applied Biosystems). 

Of the 10 loci screened, one did not amplify a product 
and one was monomorphic. The remaining eight loci 
were used to genotype and assess genetic diversity in 48 
individuals from two natural populations (Table 1).The 
Donnybrook population is situated in open Jarrah forest 
170 km south of Perth, while Lesueur is 210 km north 
of Perth in Kwongan heath. Both populations contained 
several hundred plants. Linkage disequilibrium among 
pairs of loci and deviations from Hardy-Weinberg 
equilibrium (HWE) were assessed using exact tests in 
Genepop 4.2 (Rousset 2008), applying the Bonferroni 
correction. Null allele frequencies were estimated for 
each locus using INEst 1.1 (Chybicki and Burczyk 2009) 
and genetic diversity parameters were estimated with 
Genepop 4.2. 

Results and discussion 

Evidence for null alleles and deviations from HWE was 
found for Ahum8 in both populations and for AhumlS 
at Mt Lesueur (Table 1). As no deviations from HWE 
were found at other loci and the species is dioecious, 
the observed deviations were most likely caused by null 
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Table 1. Primer sequences and characteristics of ten microsatellite loci isolated from Aflocasuanna humilis. See main text for details of PCR programs. Allele size range and number of 
alleles (N A ) are based on 48 individuals from two natural populations. Observed and expected heterozygosity (tf, H) and null allele frequencies are based on 24 individuals within each 

natural population. One locus was monomorphic and one did not amplify. 


Locus 

Primer sequence (5'-3')a 

GenBank 

accession 

Repeat 

motif 

MgCI 2 

(mM) 

PCR 

progr 

Size range 
(bp) 

N * 

Donnybrook (n = 24) 
(33°34'38.3"S, 115"51'9.2”E) 

Mt. Lesueur (n = 24) 
(30°04'03.3"S, 115” 07’40.5" E) 



Null 

allele 

H 

o 

H 

€ 

Null 

allele 

Ahum6 

F: TAAGCACGC AAAGAAACGCT 

KM374578 

(AT), 

2.5 

2 b 

104-128 

10 

0.826 

0.820 

0.034 

0.500 

0.705 

0.105 


a GGTCTTCAT'TTCAATCTCTTGGG 













Ahum8 

F: CCACTGTTA^GGCCTACAGGA 

KM374579 

(AG)„ 

2.5 

1 b 

94-122 

12 

0.522* 

0.853 

0.141* 

0.542* 

0.893 

0.144* 


fc ACACAGTCAAG ACGGCTTTC 













Ahum9 

F:TGATGAATGTAATGAAGGCTGC 

KM374580 

(GA), 

2.5 

1 

97-131 

14 

0.500 

0.553 

0.044 

0.792 

0.650 

0.024 

ft TGGATTTAAGCTACCTGATGGA 













AhumIS 

F: CCAACCAAT CTACAGAACAGCA 

KM374581 

(AT), 0 

2.5 

1 b 

186-236 

17 

0.591 

0.827 

0.096 

0.208* 

0.920 

0.272* 


RiTCTAACGTTG (GTCAAGCAATCA 













Ahumlb 

Fr GCTGTCAGC TCTGGGTAAATTC 

KM374582 

(CT), 

2.5 

1 b 

159-191 

12 

0.833 

0.789 

0.028 

0.875 

0.815 

0.025 


a TGTACACGTTGG AGTTGTGGA 













Ahum18 

F: AAGGAAATC i-CCCATAATCGAA 

KM374583 

(ATT) 20 

2.25 

1 b 

136-199 

19 

0.956 

0.896 

0.019 

0.956 

0.947 

0.019 


a AGCAAATTTGTCCAACACCC 













Ahum22 

F: CAAGGCTTC AACCAAAGATCA 

KM374584 

(AT), 

2.5 

l b 

129-141 

6 

0.667 

0.768 

0.064 

0.667 

0.576 

0.043 


a GCTGTACATGTTGGCAAGTCA 













Ahum24 

F: T1TTG AGCGTATGCAATGGG 

KM374585 

(GT), 0 

2.5 

1 b 

110-122 

7 

0.783 

0.808 

0.042 

0.583 

0.662 

0.076 


R: TCACCATTT/.AATTGAGGGTCT 













Ahum 11 

F:TTTCTTTGC7 r CGCTTTCTATGTTG 

KT224416 

(AC), 

2.5 

l b 

90 

0 








R: ACGTCAATTGATCCAAAG AAA 













Ahum14 

F: TCCCGGGATTGTAAGTGTTG 

KT224417 

(TTCTGA) 9 

2.5 

l b 

191 

1 








aTATCTCCTGCGAGACGTCTG 














a Forward primers include an Ml 3 universal sequence appended to the 5'end (5 - CACGACGTTGTAAAACGAC) 
b 0.5M betaine included in PCR 

* Significant deviation from Hardy-Weinberg equilibrium or significant estimate of null allele frequency (P < 0.05) 


Microsatellite markers for the dwarf sheoak 














































































Llorens et al. 


alleles rather than inbreeding. There was no evidence 
for linkage disequilibrium between any loci in either 
population. Among the 48 screened individuals, the 
number of alleles per locus ranged from six (Ahum22) 
to 19 (Ahum18), with a mean of 12.1. We detected 
moderate to high genetic diversity, with expected 
heterozygosity ranging from 0.553 to 0.947 (mean 
0.789) and observed heterozygosity from 0.208 to 0.956 
(mean 0.710). 

These microsatellite markers are the first to be 
developed for a Western Australian Allocasuarina and 
have doubled the number of markers available for 
the genus. Reliable polymorphic microsatellites have 
proved difficult to isolate in Allocasuarina (Lamont 
2010; Lamont et al. 2012; L. Broadhust pers. comm; 
this study). Six microsatellite markers have previously 
been developed for A. verticillata (Lam.) LASJohnson 
(Broadhurst 2011) and two for A . emuina LASJohnson 
(Lamont et al. 2012). Consistent with the relatively 
recent radiation of the genus (Steane et al. 2003), the 
markers developed by Broadhurst (2011) and Lamont 
et al. (2012) had high rates of transferability among 
Allocasuarina species both within and across sections 
(Lamont 2010). The eight micro satellite markers 
developed here for A. humilis are therefore expected to 
be transferable to other Allocasuarina species for which 
assessment of genetic characteristics may be required 
for conservation or other purposes. These markers are 
being used to evaluate the patterns of genetic diversity, 
divergence and structure among A. humilis populations 
in the SWAFR in order to inform appropriate seed 
sourcing strategies for restoration of disturbed or 
degraded land. 
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